Non-perturbative approach to backscattering off a dynamical impurity in ID Fermi 

systems 
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We investigate the problem of backscattering off a time-dependent impurity in a one-dimensional 
electron gas. By combining the Schwinger-Keldysh method with an adiabatic approximation in 
order to deal with the corresponding out of equilibrium Dirac equation, we compute the total energy 
density (TED) of the system. We show how the free fermion TED is distorted by the backscattering 
amplitude and the geometry of the impurity. 
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The physics of tunneling through static barriers is a 
topic of great importance in the subject of correlated 
quasi-one dimensional electron transport. A considerable 
number of well-established results derived over the last 
few years constitute by now the standard knowledge in 
the fieldi. Relevant applications of the theoretical find- 
ings comprise the behavior of strongly anisotropic physi- 
cal systems such as organic conductors 2 :, charge transfer 
salts^, quantum wires^, carbon nanotubes^ and quantum 
Hall junctions-. 

More recently, some attention has been paid to the 
problem of electron transport through time-dependent lo- 
calized perturbations in one dimensional correlated mat- 
ter. The interest has been to a large extent driven by the 
possibility of pumping charge and spin into a conduc- 
tor (or semi-conductor) by means of an induced effective 
time-dependent potential^. Actually, it is not difficult 
to devise examples of relevant experimental setups which 
would provide the realistic arena for time-dependent one- 
body interactions. A laser beam applied to a carbon nan- 
otube is likely to produce a periodic deformation of the 
lattice structure, playing the role of an effective time- 
dependent impurity; a similar effect should be expected 
in the localization of optical phonons in topological de- 
fects. It is worth noting, furthermore, that the clean- 
est system for experimental investigation would be given, 
probably, by a quantum Hall bar with a time-dependent 
gate producing a harmonic backscattering between edge 
states. 

Despite the possible technological advances that the 
control of charge and spin currents could bring up, we are 
faced with the fundamental question of what we could 
learn about electron correlations using time-dependent 
potentials as out-of-equilibrium probes. To answer this 
question, at least partially, more theoretical work is 
needed to understand the detailed dynamics induced by 
this type of perturbation. 

A potentially interesting observable that characterizes 



a tunneling process is the energy resolved current j(u>). 
In Ref. HI, the relevance of this quantity in the upper 
region of the spectrum (i.e., above the Fermi energy, 
lu > Ep) was emphasized, if one is interested to get in- 
formation about correlations in the leads. To leading or- 
der in the tunneling amplitude, the energy resolved cur- 
rent j{ui) is related to a simpler observable, namely, the 
electron energy distribution function n(uj), also referred 
to as the total energy density (TED) in the literature 9 
(for a precise mathematical definition see Eq. (JHJ be- 
low). In general grounds, n(u)) gives information about 
the perturbation of the ground state in the leads due 
to tunneling processes. For a non-correlated material, 
for instance, n(ui) should vanish above the Fermi sur- 
face. Any population of the spectrum above the Fermi 
energy is originated from a combined effect of correla- 
tions and multi-particle tunneling, due, on its turn, to 
out of equilibrium processes^. In Ref. 10, n(oj) was eval- 
uated in a model of correlated one-dimensional fermions 
with a time-dependent impurity coupled to the elec- 
tron density through a forward-scattering coupling. In 
Ref. [ll| this model was analyzed by means of functional 
bosonizationi^, focusing, in particular, on the transients 
produced by turning on the oscillatory impurity strength. 

It is important to notice, however, that backscattering 
effects are expected to be relevant as a rule, in all but 
rather exceptional experimental settings 13 . Recently 14 , 
effects of backscattering in a time-dependent ultralocal- 
ized impurity were studied perturbatively, suggesting a 
striking enhancement of the total current for special val- 
ues of the Luttinger parameters. Although very inter- 
esting by itself, it is necessary to take some care with 
perturbative calculations of tunneling processes, since 
multiple tunneling events may be missed in the series 
expansions, specially when dealing with finite barriers. 
The problem of backscattering by dynamical impurities 
is usually a very difficult one. Some models, like the 
spinless Luttinger model with a delta-like impurity, can 
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be solved exactly for the specific value of the Luttinger 
parameter K = f/2 1 ^. However, for general strengths of 
the electron-electron interaction and finite-ranged impu- 
rities, there are no available closed analytical solutions 
(even in the free case). Therefore, it is of fundamen- 
tal importance to address non-pcrturbative strategies in 
such a context. 

In this work we develop a non-perturbative calculation 
for a somewhat simple model, chosen to illustrate our 
method. We are interested in studying the effects that 
the backscattering off an extended dynamical impurity 
of width a, oscillating with frequency f2, will have on 
the spectrum of a one-dimensional fcrmion gas. For this 
purpose we have implemented an adiabatic approxima- 
tion which allows one to compute the TED in a straight- 
forward way. The time scale of the barrier oscillations 
is given by I/O, while on the other hand the traversal 
time for tunneling is given by a/\vl^, where v is the ve- 
locity of fermions inside the barrier. Therefore, if the 
traversal time is much smaller than the oscillation time 
(a/\v\ <C 1/S1), the tunneling can be considered as tak- 
ing place through an essentially static barrier. Thus, the 
adiabatic approximation just consists in the calculation 
of the spectrum or, in general, any correlation function, 
in the limit f2a/|i;| — > 0, neglecting subleading correc- 
tions, of order (Via/\v\) 2 . Interestingly, the adiabatic 
regime was recognized as the relevant one in the con- 
text of charge quantization in pumping processes^. It 
is correct to state that in the adiabatic limit it is possi- 
ble to evaluate in an exact way any fermionic correlation 
function without relying neither on a perturbative ex- 
pansion in the backscattering coupling constant nor on 
the finite range of the dynamical barrier. This means, in 
principle, that it is possible to capture multiple tunneling 
processes and bound states, which are absent in the case 
of ultra-local potentials (<5-like potentials) and small val- 
ues of the coupling constants. In appendix (QQ) we have 
applied the method to the exactly solvable case of pure 
forward scattering to explicitly show how the adiabatic 
approximation works. 

As the computational starting point, let us consider 
the following Hamiltonian, which describes the interac- 
tion of spinless fermions with an external effective time- 
dependent potential V(x,t), responsible for backscatter- 
ing transitions between right and left movers: 



H — H + H imp , 



(1) 



where 



H Q =ihv F / dx Otp R d x i/} R - V'l^xV'i) > ( 2 ) 



and 

H imp =g b J dx {i) R i)h + V'iV'il) V(x,t) . (3) 

Above, gb is the coupling constant associated to the back- 
ward scattering of electrons caused by the presence of a 



time-dependent harmonic barrier. It is also possible to 
consider a more general Hamiltonian with forward scat- 
tering couplings in addition to the backward coupling 
considered here, however, this does not lead to any ad- 
ditional interesting physical effect. Our results can be in 
fact extended without much difficulty to the case where 
a forward scattering coupling is taken into account. 

Although we have verified that both the method and 
the general results are independent of the explicit details 
of V{x, t), we will use, just to fix ideas, a square potential 
profile, 

V{x, t) = (0(x + a/2) - @(x - a/2)) cos(fti) , (4) 

where a is the width of the square potential and f2 the 
oscillation frequency. 

We are particularly interested in obtaining the TED 
for the above model. We recall that in the Wigner repre- 
sentation the TED can be written in terms of the fermion 
correlation function as 

/oc 
dr e ZbJT G+- (r = 0, R, T, r) , (5) 
-oc 

where we have introduced the closed time path 
formalism 18 in which fermion propagators are time- 
ordered along the usual Schwinger-Keldysh time contour 
C: 



G c = 



G++ Ga 



G — l G — 



(6) 



The subscripts + and — refer to fields defined in the 
upper and lower branches of C, respectively, correspond- 
ing to the forward (+) and backward (— ) time evolution. 
Above, r, r and R, T are the spatial and temporal relative 
and center of mass coordinates, respectively. 

The static case, il = 0, is exactly solvable. It corre- 
sponds to take the limit in Eq. JSJ) such that, 



lim n(u>, R - 



0,T) 



(7) 



where we choose to calculate n sta tic(w) at the center of 
the barrier R — 0. 

Physically, the above limit means that the TED is as- 
sumed to be probed in very short time scales, involving 
an external process with a relaxation time To <C I/f2. In 
this regime, the TED behaves as effectively static and is 
given by Eq. (0. 

On the other hand, if T > l/ft, the TED given by 
Eq. 10 is a very rapid oscillating function of time and 
it is not the relevant quantity to evaluate. An arbitary 
observable in this regime should be computed by means 
of an averaging procedure. The meaningful physical ob- 
servable in our case is, therefore, the averaged TED 



n(uj,R = 0,T)dT , 



(8) 
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It is important to notice that this averaged quantity is 
purely dynamical and is in general disconnected from the 
"instantaneous" or static regime. Although we can take 
the limit — ■> in Eq. (JSJ, this does not correspond 
to the static limit since, the definition of the averaged 
TED implies that it should be probed over an infinite 
time interval. Mathematically, it is simple to understand 
from Eqs. (J7J and iJSJ that 



lim n(oj) ^ n sta ticM 



(9) 



once we cannot interchange the limit with the integral 
operator. 

Since we arc interested in the dynamical aspects of the 
problem we concentrate ourselves in the calculation of 
the averaged TED given by Eq. For this purpose, we 
have calculated the Green's function and TED, Eq. JSJ), 
in the adiabatic approximation, which, for the potential 
given by (0J amounts to substituting gtcos(flt) — > gt,. 
Then, we have restored the temporal dependence and 
evaluated the averaged TED using Eq. ||SJ). The validity 
of this procedure (i. e. of restoring the time dependence 
in n(to) after integrating the Green's function), relies in 
the extra condition oj 3> £7, which allows a practical com- 
putation of the TED, and preserves non-trivial physics, 
holding beyond perturbation theory. The relaxation of 
this condition forces us to restore the time dependence 
at the Green's function level, turning the calculation con- 
siderable more difficult. Therefore, our results are useful 
in the frequency range u> 3> f2 3> 1/Tq. The first inequal- 
ity is imposed to simplify calculations, while the second 
one is implicit in the definition of the averaged TED, Eq. 
©. 

Our main results are depicted in Figs. n an d|3 where 
n(uj) is shown for fixed gb, with variable a (Fig^l, and 
for fixed a, with variable gb (Fig|5J). Observe that at this 
level of approximation the curves are independent of f2. 
It is worth noting, however, that the range of validity of 
our predictions does depend on the external frequency 
due to the condition u> 3> f2. Also, in both cases, in con- 
trast to the situation in which only forward scattering is 
taken into accounliSiii, no gain or loss of energy-quanta 
of nfl takes place (see appendix [XJ . This is in fact a 
consequence of the oj ^> regime considered here. How- 
ever, one may clearly observe that backscattering breaks 
the uniformity of the Fermi sea, characteristic of the free 
electron gas. The TED n(oj) has its maximum peak for 
oj < —\gb\- For oj > —\gb\, the TED shows a pronounced 
decay, if compared to the free case behavior. When \agb\ 
grows, the peak tends to be situated at w ~ — \gb\, and 
the valley observed for oj > — \gb\ is drastically lowered, 
indicating an important depletion in the population near 
the Fermi surface (in this calculation we have tuned the 
Fermi energy to oj — 0). Concerning the main peak, it 
is interesting to note that the maximum value of the av- 
eraged TED obeys the power law max\n(oj)] « y/\ag b \, 
for \agb\ sufficiently large. Although the present model 
does not contain correlations, it is clear from these results 
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FIG. 1: Averaged energy density n(oj) as a function of the 
frequency oj for fixed backscattering coupling and variable 
potential range a. oj and gb are measured in units of the Fermi 
energy Ef = vfIzf while the potential range a is measured 
in units of the inverse of the Fermi momentum l/kp. 
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FIG. 2: Averaged energy density n{ui) as a function of the 
frequency oj for fixed potential range a and variable backscat- 
tering coupling g b . 



that the presence of backscattering off the dynamical im- 
purity will modify correlations in a relevant way, mainly 
for large values of \agb\, when the structure of the elec- 
tronic density near the Fermi surface is strongly affected. 
Let us emphasize again that the peak structure we found 
(for the region lo ^ O ^> 1/Tq) is not of static origin. 
Indeed, the static TED can be easily computed (Eq. Q) 
and it presents an oscillating structure with maxima of 
equal height. 
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We now sketch some technical details concerning the 
adiabatic approximation. The main idea of this approxi- 
mation, as stated above, is that in a regime of sufficiently 
slow barrier oscillations, Oa/|u| <C 1, the spectrum of the 
system can be considered essentially static. In this case, 
we can safely consider (formally) <7&cos(f2t) ~ g b . The 
Dirac equation corresponding to the Hamiltonian Q can 
be written down as 



id t ^ = (H a + H imp )^. 



(10) 



We seek, then, for stationary solutions of energy E. 
Solving the resulting pair of coupled differential equa- 
tions for the right and left components of '5, we get for 
the right field (to save notation we drop the index R): 



1 



-ia(E+k)/2 



Z7r E+k 



ikx 



gb e 



E + k 

in the region —a/2 < x < a/2, where k 2 = E 2 



(11) 



5 b 2 (fr 



now on, we take h = l,vp = 1). We obtain a similar 
result for the left moving fermion. Of course, outside 
the barrier, where V(x,t) — 0, one has just plane wave 
solutions for both chiral components. Inside the barrier, 
due to the presence of backscattering, the solutions are 
plane waves or exponential decays, depending on the sign 
of k 2 . This, in turn, determines the energy regions in 
which one or the other kind of wave function is defined. 

In terms of the wave function obtained above, the right 
moving component of the Green's function reads: 



G + -{R,r,T,r) =i [ dE 

J — OO 



** E (R - r/2, T - t/2)V e (R + r/2, T + r/2) . (12) 

Substituting Ijllfl into l|12l) . using Eq. JSJ and explicitly 
evaluating the integrals, we obtain for the TED at the 
center of the barrier the following expression: 



, j = (e(-^)-9(-^-|fffr|))(ff& + ^) + Q(-uj-\g b \) (g b +u+^/u 2 -g 2 ) 2 + a 

cosh(y/ g 2 - uj 2 a) g b + lu (uj + ^lu 2 - g 2 ) 2 + g 2 + 2g b {u + ^uj 2 - g 2 ) cosh^^ _ w ^ a ) 

. (131 



The first term in the right hand side of Ijl3(l comes from 
the contribution of non-plane wave solutions (\E\ < \g b \), 
while the second term comes from the energy region 
(\E\ > \g b \) where the spectrum is composed essentially 
by plane waves. It is interesting to note that no pertur- 
bation theory in the backscattering amplitude g b is able 
to capture the physics of the first term, since perturba- 
tion theory only deals with small perturbation on plane 
waves. 

Finally, restoring the temporal dependence through 
the formal substitution g b — > g b cos(flT), and perform- 
ing the average over a period (Eq. (JSJ), we obtain the 
plots shown in Figs. 1 and 2, discussed above. 

To summarize, we have studied the effects of backscat- 
tering in a one dimensional fermion gas with an oscil- 
latory square barrier. We have calculated the TED in 
the adiabatic approximation. This method allowed us to 
consider finite barriers in a non-perturbative regime. In 
this way we were able to obtain the dependence of TED 
with the coupling constant g b and with the geometry of 
the barrier, given essentially by its width a. We have 
found that the structure of the Fermi sea may be strongly 
modified by the presence of the dynamical barrier, pro- 
ducing a peak structure in the TED at energies around 
the typical backscattering energy g b . Also, the probabil- 
ity of finding electrons near the Fermi surfaces may be 
drastically suppressed when the parameter ag b becomes 
large enough. We also found an interesting power law 
dependence of the maximum of the TED peak, given by 



max\n(u>)] y/ 1 ag b | . This structure opens the interest- 
ing possibility of the experimental determination of mi- 
croscopic quantities like the backscattering strength g b 
and the effective width of the potential a. Of course, the 
spectrum modifications are expected to affect the elec- 
tron correlations in the wire not trivially. We are cur- 
rently analyzing the combined effect of electron-electron 
interactions in the wire with backscattering due to strong 
and extended dynamical barriers. 
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APPENDIX A: FORWARD SCATTERING 

In this section we consider an exactly solvable problem, 
in order to ilustrate how the adiabatic approximation 
works. 
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Let us consider the Hamiltonian 



(^) 2 "e(-u;-Q(n-i)) 



(A4) 



H = H Q + H i 



imp ) 



(Al) 



where 



Hq = i h vp I dx 



(fad^R-lid^i) , (A2) 



and 



#imp = 3/ J dx (lp R 1pR 



+ rM v( X ,t) 



(A3) 



Above, <7/ is the coupling constant associated to the for- 
ward scattering of electrons caused by the presence of a 
time-dependent harmonic barrier. 

In this case, the Dirac equation (|10fl is easily solved, 
since the right and left components of the spinor decou- 
ple. Then, the right mover component is given by a ex- 
pression of the type ipR ~ e%kx while the left one reads 
tpL ~ e~~ lkx . Now, it is a simple matter to calculate the 
Green's function (Eq. (|T2l ') and the averaged TED ( Eq. 
iJHJ) is given in the adiabatic approximation by, 



oo 2n 



*»m = EEt- 1 ) 4 - 



(2n) 



n=0 ; 



n\) 2 {2n-i)W. 



where we have chosen h = 1 and = 1. 

However, in this case, it is not necessary to invoke 
any type of approximation, since this problem is ex- 
actly solvable 10,11,12 . The reason behind this, is that 
this decoupling property between right and left movers 
is present in the full quantum problem and not merely in 
the static Dirac equation. The exact result for TED in 
the forward scattering case is given by 



n[oj 



) = E St-*) 1 



(2n)! 



n=0 i=0 



(n\) 2 (2n-i)\i\ 



gf sin[fia/4] 



2n 



e{-uj - ft(n - i)) . (A5) 



Comparing Eqs. (|A4I) and (| A5|> . we clearly verify that 
the adiabatic approximation gives the correct result at 
leading order in 2a <§; \ (here, we have recovered the 
original units). 
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